New topological surface state in layered topological insulators: unoccupied Dirac cone 
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The unoccupied states in topological insulators Bi2Se3, PbSb2Te4, and Pb2Bi2Te2S3 are studied 
by the density functional theory methods. It is shown that a surface state with linear dispersion 
emerges in the inverted conduction band energy gap at the center of the surface Brillouin zone on 
the (0001) surface of these insulators. The alternative expression of Z2 invariant allowed us to show 
that a necessary condition for the existence of the second T Dirac cone is the presence of local gaps 
at the time reversal invariant momentum points of the bulk spectrum and change of parity in one 
of these points. 
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PACS numbers: 73.20.-r, 71.70.Ej, 72.25.Dc 

Three-dimensional topological insulators (TIs) are 
characterized by presence of the topological surface state 
(SS) emergence in the principal band gap [ljj. This 
spin-polarized linearly dispersing surface state forming 
Dirac cone arises from a symmetry inversion of the 
bulk bands at band gap edges owing to the strong 
spin-orbit interaction (SOI). Time reversal symmetry 
protects topological SS from backscattering in presence 
of weak perturbation j2)-j3] . This causes new possibility 
for practical applications particularly in realization of 
dissipationless spin transport which can be used in new 
spintronic devises. 

A number of materials that hold non-trivial 
spin-polarized Dirac state has been intensively studied. 
Among them the Bi and Sb chalcogcnidcs arc most 
investigated TIs at present. Their band structure 
is widely studied 0, i-ES] and spin texture of the 
topological SS was observed both indirectly by use of 
circular dichroism II ill and directly by spin-resolved 
experiments [B], [T3, Il3^ . However, besides the Dirac cone 
various types of SS states take place in TIs. Some of them 
can result in new physical phenomena. For example, 
angle-resolved photoemission spectroscopy (ARPES) 
measurements as well as time- and angle-resolved 
photoemission spectroscopy experiments |l4[ led to 
discovery of parabolic spin-split surface states in the 
principal energy gap just below the conduction band and 
M-shaped states in the local valence band gap 



Here we study the conduction band (CB) energy 
gap surface states in the layered topological insulators 
Bi 2 Se 3 , PbSb 2 Te4, and Pb 2 Bi 2 Te 2 S3 by using density 
function theory (DFT) methods. The choice of these 
TIs is motivated by the fact they have different crystal 
structure, composed of quintuple (Bi 2 Se3), septuple 
(PbSb 2 Te 4 ), and nonuple (Pb 2 Bi 2 Te 2 S 3 ) layer blocks. 



Bi 2 Se3 is the extensively studied TI while PbSb 2 Te4 and 
Pb 2 Bi 2 Te 2 S 3 were recently predicted as TIs [U,!!!- We 
show that except the Dirac cone in the principal gap and 
trivial unoccupied surface states a massless spin polarized 
surface state arises in the CB local gap at the T. This 
state induced by the SOI inverted local gap has spin 
hclicity similar to that of the Dirac state in the principal 
energy gap. 

For electronic band calculations we employ two 
different computer codes that are based on DFT. The 
first one is the Vienna Ab Initio Simulation Package 
(VASP) [25], HH- We used the generalized gradient 
approximation (GGA) [2?| to the exchange correlation 
potential and the projector augmented wave (PAW) 0, 
|29( basis sets to solve the resulting Kohn-Sham system. 
The second approach used for electronic structure 
calculations is the full-potential linearized augmented 
plane-wave (FLAPW) method as implemented in the 
FLEUR code [13| with PBE for the exchange-correlation 
potential. The FLAPW basis has been extended 
by conventional local orbitals to treat quite shallow 
semi-core (i-states. Both methods contained scalar 
relativistic corrections and spin-orbit coupling was taken 
into account by the second variation method [3l|]. To 
simulate the (0001) surfaces of the TIs we use a slabs 
composed of 6 quintuple layers (QLs) and 6 septuple 
layers (SLs) for Bi 2 Se3 and PbSb 2 Te4, respectively, and 
5 nonuple layers (NLs) for Pb 2 Bi 2 Te 2 S3. 

The calculated unoccupied electronic structure spectra 
in materials under interest is shown in Fig. [TJ Shaded 
regions depict projected bulk bands onto the (0001) 
plane, while red dots indicate a weight of the states in the 
outermost structural block (QL, SL, NL). As easy to see 
in addition to the Dirac cone in the principal gap there 
exist different trivial surface states at the energies shown 
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FIG. 1: Fig. 1 Surface electronic structure of Bi2Se3 (a), PbSb2Te4 (b), and Pb2Bi2Te2S3 (c) above _Ep as calculated by VASP. 
Size of red circles correspond to the weight of the states in the outermost block. The projected bulk bands is shown in green. 
Dashed frames mark region of the unoccupied Dirac-like surface state. 



in Fig. [T] and at higher energies in the local CB gaps in 
all surfaces of interest. 

At the same time in the conduction band energy gap 
there appears a Dirac-like surface state at the T point 
(framed in Figs.QJa-c) by dashed line) for all compounds 
we study. In the case of Bi2Se3 this state arises at ~ 1.5 
cV in gap of 70 meV formed by projection of the second 
and third bulk conduction bands. For PbSb2Te4 and 
Pb2Bi2Te2S3 the Dirac-like surface state takes place at 
~ 0.8 eV and ~ 1.8 eV respectively. These states have 
the linear dispersion in the local gap and entering to the 
bulk states projection they transform into resonances and 
spread along the edges of the projection. A close-up 
of these states is shown in Fig. [5J Note that the 
Dirac-like surface state in the conduction band is well 
reproduced by both VASP and FLEUR calculations with 
some differences in the projected bulk bands (Figs. [2] (c) 
and (d)). 

In all considered TIs the unoccupied Dirac-like SS 
is isotropic with respect to fen, i.e. the constant 
energy contours below and above the degeneracy point 
have an ideal circular shape as schematically shown 
in Fig. [3] The calculated orientation of the electron 
spin within the surface state demonstrates in-plane spin 
polarization (out-of-plane component is negligibly small) 
with positive (clockwise) spin helicity in the upper part 
and negative helicity in the lower part of the Dirac cone, 
i. e. the same helicity as in the case of the topological SS 
located in the principal energy gap. Besides, owing to the 
SOI entanglement of the spin and orbital momenta the 
spin polarization in the CB cone is reduced and does not 
exceeds 60% in Bi2Se3 and Pb2Bi2Tc2S3 and is slightly 



smaller in PbSb2Te4 (~ 50%) that is comparable with the 
spin polarization in the conventional Dirac cone fl3l l32j . 

In order to reveal the origin of formation of the second 
Dirac state in the conduction band of TIs consider first 
the bulk band structure of Bi2Se3. Note that a direct 
consequence of the inversion of the Bi and Se states in the 
principal gap at the T point in Bi2Sc3 is formation of the 
Dirac cone in the principal energy gap. Let us consider 
the bulk conduction band of Bi2Se3 in the energy range 
of interest. Without spin-orbit coupling included the 
second and third conduction bands are degenerate along 
the T — Z direction (Fig. Eta)). Spin-orbit interaction lifts 
this degeneracy and opens a gap between these bands at 
the T point (Fig. BJb)). Thus, the gap supporting the 
new massless state in the conduction band results from 
the SOI. Both the upper and lower bands are mostly 
composed of the bismuth states, however, the Sc p z states 
of the central atomic layer of QL contribute to these 
bands too. As one can see in Fig. E(c), the weight of 
Se states in the lower band goes to zero approaching the 
r point. At the same time Se states in the upper band 
have maximal weight in the vicinity of T. Thus, this 
SOI- induced local gap at ~ 1.5 cV between the second 
and third conduction bands has inverted Se states. Like 
to the principal P gap where the states of the outer Se 
atoms of QL are inverted the states of the inner Se atom 
have a similar SOI-induced inversion in the CB gap. This 
fact points out that the CB Dirac cone is the unoccupied 
topological SS. 

In the bulk band spectra of PbSb2Te4 the similar 
SOI-induced inversion can be detected for states of the 
inner Tc atoms of SL in the vicinity of the Z point 
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FIG. 2: Fig. 2 A magnified view of dashed frames marked 
in the Fig. [1] with CB cone in the center: (a) Bi2Se3; 
(b) PbSb 2 Te 4 ; Pb 2 Bi 2 Te2S3 as calculated by VASP (c) and 
FLEUR (d) codes. 



(Fig. SJd)). However, in this ternary compound the 
picture is more complicated and it becomes even less 
obvious for quaternary TI Pb2Bi2Te2S3. 

From this point of view it would be desirable to have 
such a simple technique to qualify unoccupied topological 
states as calculation of I2 invariant. The difficulty 
is that the unoccupied cone is appeared in the local 
gap of the bulk conduction band where along a path 
between time-reversal invariant momenta (TRIMs) in the 
Brillouin zone the gap in the spectrum of the bulk states 
closes and opens. 

Analysis of the relation (-1)"° = (-1) 2P3 (H 
of the index Vq 0, HH in the theory of topological 
insulators with the winding number 2P3 or with the 
non-Abelian Chern-Simons term in the fc-space [35j using 
the continuous approach in the theory of topological 
invariants does not solve t he p roblem. Thus we address 



the approach of Refs. [36H38I ] enabling us to determine 



the mechanism of reduction of the contribution made by 
bulk states to the topological invariant. 

The origin of reduction of the contribution of the 
bulk states can be found by comparing the continuous 



FIG. 3: Fig. 3 Schematic view of the spin structure of the CB 
cone. Red and blue arrows indicate spin direction for upper 
and lower part of the cone, respectively. Scales for constant 
energy contour and 2D Brillouin zone are different (typical 
maximal radius of the contour in materials under study is 
~ 0.03A" 1 ; f — K direction is of ~ lA -1 ). 
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FIG. 4: Fig. 4 Bulk band structure of Bi2Se3 without (a) 
and with (b) SOI included; (c) magnified view of blue dashed 
frame marked in the Fig. |4jb) with composition of the second 
and third conduction bands in the vicinity of P; (d) the same 
for PbSb2Te4 in the vicinity of Z. 



and lattice versions of the Z2 invariant. It is shown in 
Ref. [3(1 | that the Z2 invariant in the continuous case is 
expressed as 
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where B~ = [—ir, 7r] ® [— 7r, 0] is half of the Brillouin zone, 
A = Trip^dip and F = dA are the Berry gauge potential 
and the associated field strength respectively and tf>(k) 
is the 2M(k)-dimensional grou nd state multiplet. The 
lattice analog of Eq. (TTJ) is [37| 
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since £ fceB - F ( k ) = J2kedB- A ( k ) + 2 ^J2keB- n ( k )- 
n(k) in this equation are integers and n(fc)mod2 6 Z 2 
due to the residual U(l) invariance [37| ■ From Eq. ^) we 
can conclude that the reason for cancelation of the bulk 
continuum states is the compactness of the lattice gauge 
theory. In all TIs under study at energy of CB Dirac 
state there are gaps in T, Fj, Lj, and Z TRIMs in the 
case of rhombohedral Bi2Se3 and PbSb2Te4 or in V, Mi, 
Li, and A TRIMs of hexagonal Pb2Bi2Te2S3 which allow 
us to implement the parity analysis 34[ to calculate the 
topological number Uq for the CB cone in these systems 
by fixing of an isoenergy level in the local gap, supporting 
the unoccupied Dirac cone. 

The results for PbSb2Te4 show that non-trivial value 
Vq = 1 of the topological invariant for the unoccupied 
cone is provided by the change of the parity at the 
Z point (where we have revealed the inversion of Te^n 
states, see Fig. Hfd)). In the case of Pb2Bi2Te2S3 the 
change of parity in the A point in the gap between 4-th 
and 5-th bulk conduction bands determines topologically 
non-trivial isq and the appearance of the unoccupied 
Dirac cone in the surface electronic spectrum. Thus 
inversion of the local gap edges in the conduction band 
along with the presence of non-inverted gaps in other 
TRIMs of the bulk Brillouin zone at the same energy is 
responsible for the emergence of single Dirac-like helical 
spin surface state in the topological insulators. 

In summary, we have examined the unoccupied 
electronic spectra of the layered topological insulators 
with different composition and crystal structure. We 
revealed that except for trivial conduction band surface 
states a massless topologically-protected helical spin 
state exists in the narrow local conduction band gap 
at the r point in all materials under study. This state 
arises owing to inversion of the certain bulk conduction 
bands at the time-reversal invariance momentum points 
r, Z, and A of Bi 2 Se 3 , PbSb 2 Te 4 , and Pb 2 Bi 2 Te2S3, 
respectively. The observed state has the same spin 
helicity as in the conventional Dirac cone in the principal 
energy gap. The revealed CB topological states provide 
a pathway to the measurements on excitations between 
two Dirac spin polarized states at the surface of layered 
topological insulators. 
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